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Abstract 

We present the exact solution to a one-dimensional multicomponent quan- 
tum lattice model interacting by an exchange operator which falls off as the 
inverse-sinh-square of the distance. This interaction contains a variable range 
as a parameter, and can thus interpolate between the known solutions for the 
nearest-neighbor chain, and the inverse-square chain. The energy, susceptibil- 
ity, charge stiffness and the dispersion relations for low-lying excitations are 
explicitly calculated for the absolute ground state, as a function of both the 

range of the interaction and the number of species of fermions. 
03.65.Db 
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I. THE LATTICE MODEL FROM THE EXCHANGE MODEL 



In a recent Physical Review Letter to be called S2, Sutherland and Shastry intro- 
duced a one-dimensional, continuum, multi-component quantum many-body system with a 
hyperbolic exchange interaction, governed by the following Hamiltonian: 

Hc ° nt = ~ 2 ? dx] + £ sinh 2 ^-^)- (1) 
We have set H = 1, the mass to be unity, Pjk is the exchange operator which permutes 
particles j and k, and we always have in mind the thermodynamic limit, when the number 
of particles iV and the box size A both become infinite, while the density N/A = d remains 
finite. In S2, scattering in this system was shown to be non-diffractive, and so the system 
could be solved exactly (in the thermodynamic limit, with finite size corrections being expo- 
nentially small in A) by the asymptotic Bethe ansatz |J. The system is characterized by the 
number of different types b(f) of bosons (fermions) in the system; we write this as a B b F^ 
system. Thus B(F) represents a single component system of identical bosons(fermions) 
interacting with an interaction strength A(A =F 1). 

This system has a classical limit, obtained by taking A — > oo, and in the ground state, 
the particles freeze into a lattice with lattice spacing A/N = 1/d, so that Xj — > j/d. Ex- 
panding the Hamiltonian in 1/A, we find H cont — > H e[as + 2XHi att , where H e[as represents 
the Hamiltonian for the elastic degrees of freedom — phonons and solitons — identical to 
the single component F system ||. The additional term of order A gives the Hamiltonian 
for the compositional degrees of freedom Hi att , 

Hlatt = ^ k s^\u%/dy (2) 

There is exactly one particle to a site. Thus, to order A, the elastic and compositional 
degrees of freedom separate. Coupling between the two occurs through an expansion of the 
interaction about the regular lattice, and so is of order (Sx 2 ) ~ 1/A, and down by another 
factor of 1/A. (A similar procedure is used to derive the Heisenberg antiferromagnet from 
the half-filled Hubbard model.) 



In this letter, we will study the spectrum of this lattice Hamiltonian. The Hamiltonian 
has many realizations — as does the corresponding nearest-neighbor model — such as by 
spins, a (supersymmetric) t-J model, SU(N) spins, etc. In fact, this system interpolates 
between the nearest-neighbor model ||,|] (d — > 0), and the 1/r 2 lattice (d — ► oo) 

treated in S2. Our solution here will rely heavily on Sutherland's treatment || - to be 
called SI — of the nearest-neighbor system. Two results proven in SI can easily be shown 
for the hyperbolic case: (i) the ground state of the B b F^ system is identical to the ground 
state of the BF? system, so for the ground state, we need never consider more than one type 
of boson which we usually think of as a vacancy; (ii) the Hamiltonian with Pjk — > —Pjk is 
unitarily equivalent to Hi att , hence "ferro-" and "antiferro-" cases are both included in the 
ground state of Hi att . The pure (anti)ferro- case is B b {F^). Partial results concerning the 
integrability of this model have previously appeared in the literature ||. 

II. THE SOLUTION OF THE EXCHANGE MODEL 

Because the system is non- diffract ive, the wavefunction must be given asymptotically as 

JV 

V(x\Q)~Y,nP\Q)ex9[iY,*i k Pi\> ( 3 ) 

P 3=1 

where Xi < . . . < x^, P is one of the N\ permutations of the N asymptotic momenta kj, Q 
is one of the N\ rearrangements of the particles, and \fr(P\Q) are (N\) 2 amplitudes related 
by two-body scattering. The fc's are determined by a set of eigenvalue equations, which say 
that the phase shift of a particle going around the ring must be unity mod 2n: 

exp [ikjA] J] - ex P H°o( k j ~ k m)} Xjj-! ■ ■ ■ X jtl X j>N ■ ■ ■ X j>j+1 * = (4) 

Here Xj jTn = X(kj — k m ) are two-body scattering operators identical to those for the 5- 
function gas or the nearest-neighbor lattice model; the eigenvectors and eigenvalues are given 
in SI. The two-body phase shift 0o(k) is the phase-shift for scattering from a 1/ sinh 2 (r) 
potential, and is given as 
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o (k) = Hog 



r(l + ifc/2)r(A — ifc/2) 



(5) 



_r(l — i/s/2)r(A + ifc/2) 
Using the eigenvalues of the X's from SI, one obtains coupled equations for the fc's, the 
first being 



Lk = 2nl(k) + J2 9 o( k - k )-J2 9 ( 2 ( k ~ «))■ 



(6) 



Here 9(k) is the phase-shift for a 5-function gas: 9{k) = 2 arctan(/c/2A). The I(k) are the 
quantum numbers, from log(±l)/27r. The remaining equations are identical with those of 
the nearest-neighbor model; the fc's couple directly only to the second equation, which has 
the form 



= 2nJ(a) - ]T e ( 2 ( a - k )) + 



(7) 



The remaining terms and equations are the same as for the nearest-neighbor model, and 
their form and number depend upon what types of particles we have. 



III. THE LIMIT A ->■ oo 

When we take the limit A — > oo, we must also rescale the asymptotic momenta k, defining 
a new variable x = k/2X. Likewise the asymptotic momenta a for the compositional degrees 
of freedom must be rescaled, but we use the same symbol, so a — > a/2 A. Then the phase 
shift for particle-particle scattering has an expansion in 1/ A of the form 

e (k)^2\e (x) + e l (x) + ..., (8) 

where 

6 (x) = \\x log(l + l/x 2 ) + % log[(l - ix)/{l + ix)] \ , (9a) 

»JW = bog(l+ l/l 2 ), (9b) 

^ = -\ + ih- < 9c > 

Finally, for one more redefinition, let 6(x) — > 2arctan(x). 



We now expand the first equation in 1/A, obtaining to order A — which we call the 
zeroth order, since it starts the expansion — the equation 

A X = J2%(x -x). (10) 

x 

Let Np(x)dx be the number of x's in x — > x + dx. Then the equation for p becomes an 
integral equation, and upon differentiating 

1 1 r A 



9' (x -x)p(x)dx = Lp. (11) 

Zna Z7T J -a 

Let L~ l be the inverse of L, so p = L~ 1 [l/2nd]. The normalization of p is 

rA 

1 = / p( x )dx = lL + p, (12) 

J— A 

while the classical ground state energy is 

W 1 cA i oo -r 

— ^- = -/ x 2 p(a;)da; = (a^/2)V=-y — . (13) 

4A 2 iV 2J-A HK J y 1 >H 4^sinh 2 (j/d) V ; 

Before discussing the solution to this integral equation, we go on and examine the first order 
corrections to the equation for the x's. 

IV. THE FIRST ORDER EQUATION 



We write the x's which satisfy the zeroth order Eq. ( |i0|) as xq, and then look for cor- 
rections to the x's as x = xq + 5x/2X + . . .. Let 7(2) = fe(x)p(a;), and NR(a)da = 
number of a's in da. Then the first order equation is 

L 7 = ^ + [ A 6'Ax - x)p{x)dx T 0(2(a: - a))R{a)da. (14) 

iV 27T J-A 27T 7-B 

Now this equation is linear, so we can write 7 = 70 + 7i, where 70 is the correction for the 
elastic modes, while 71 is the correction from the compositional modes. We are interested 
only in the compositional degrees of freedom in this Letter, so we write 7 for 71, which obeys 
the integral equation 
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L 7 = — - [ B 6(2(x - a))R(a)da = -g. (15) 
2% J-B 

Expanding the energy, and momentum to the same order, we find 

AE/2XN = xt 7 = -x^L-'g = -g'L^x = -g ] e = Jg = E latt /N, (16a) 
AP/dN = = -tfL^g/d = —g^L v \jd = -2ng^p = Pi att /N. (16b) 

Here, e and the derivative e 1 obey the equations 

Le = x, and Le = x 2 /2 - /i, e{±A) = 0. (17) 

Remember, all integrals over x are integrals only over the interval —A to A. These ex- 
pressions of Eq. ([TBI) are exactly the energy and momentum for the lattice Hamiltonian 
Hi a u of Eq. (fj). The classical ground state energy in Eq. C p~3|) can also be expressed as 
E/A\ 2 N = t^e/2nd + fi. 

As shown in Eq. ([?[), the only effect of the x's — and thus the fc's — on the a's is through 
the expression 



Y / 6(2(a-x)) = N [ 6(2(a-x))p(x)dx = Np(a). (18) 

x J -A 

Returning to the expressions for the energy and momentum of the lattice Hamiltonian 



Eq. flTBp, and substituting the expression of Eq. (|15D for g, we see that they can be written 



as 



Piatt/N=f p{a)R{a)da=p^BR, (19a) 

J-B 

Eiatt/N = / e(a)R(a)da = e^BR, (19b) 

J-B 



with e defined by 



e(a) = - [ A 9'(2(a - x))e(x)dx. (20) 

7T J -A 



Eliminating a between p(a) and e(a) gives the dispersion relation e(p) for a single excitation 
of the lattice Hamiltonian. This can also be expressed by the Fourier transform of the 
hopping matrix element, and so 
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This dispersion curve is shown in Fig. ([T|). 

V. THE DENSITIES p(x) AND e{x) 
As will become clearer, the solution to the hyperbolic lattice problem depends crucially 



on the two classical ground state densities p(x) and e(x), obeying Eq. (|TTD , ([T|) and (|Tj 
Unfortunately, these equations cannot be solved in closed form, even though we know the 
densities obey an infinite number of sum rules, arising from the classical limit of the Lax 
equations. The energy sum rule of Eq. ( |13D is the first of this heirarchy. 

However, we can derive what appears to be a convergent expansion for p(x) and e(x) in 
terms of the Chebyshev T and U polynomials, of the form 

1 



p(x) 



7rV 'A 2 — X 2 



(22a) 



e(x) = -2VA 2 -x 2 J2 ejU 2j {x/A). (22b) 

3=0 

In the nearest-neighbor limit of A — > 0, when 6' (x) — > — log \x\, one finds that the expan- 



sions stop at the first term, so that p(x) — > 1/ny/A 2 — x 2 ~ S(x), and e(x) — > — 2\j A 2 — x 2 ~ 
-7iA 2 5(x). This gives d' 1 -> -log[A/2]. 

At the other limit of an inverse square interaction, when A — > oo, then 6' (x) — > 7r5(x), 
so p(s) 1/2A, d 2A/n, and e -> x 2 - 2/i = x 2 - A 2 . In Fig. @ and (D we show p(x) 
and e(x) for representative values of A. 

VI. SUMMARY OF RESULTS FOR THE NEAREST-NEIGHBOR MODEL 

We make extensive use of the solution from SI for the nearest-neighbor chain, so we 
summarize the results here. We restrict ourself in this letter to the most interesting case 
of the F* system, with Nj particles of type j, and N\ > ... > Nf. Let M\ = N — N\, 



7 



M 2 = N — JVi — N 2) • • ., M/ — 1 — Nf. Then there are / — 1 coupled equations for / — 1 sets 
of roots, the jth set being Mj roots in number. For the ground state in the thermodynamic 
limit, these sets of roots distribute with densities rj(a) between the limits — bj and bj, 
normalized so that 

/bj 
rj(a)da = Mj/N = rrij, or t ] br = m. (23) 
-bj 

We have adopted a very useful notation where the / — 1 r's and m's are arranged as column 
vectors. Let K n represent an integral operator with difference kernel 

1 2r? 

K n (a) = — — — - — = 9'(na)/2n. (24) 
27r 1 + rrcr 

Finally, let £ be a column vector with components Xj = 5jiK 2 (a), b a matrix projection 
operator which imposes the integration limits rkbj on aj, and K a symmetric matrix of 
integral operators with elements K^ m = Sj tm K 2 — (5j, m -i + dj,m+i)Ki- Then the ground 
state densities satisfy the / — 1 coupled integral equations £ = r + Kbr = (1 + K)br. With 
the Hamiltonian normalized as 

H nn = -'£[1 + P j k], (25) 
j>k 

the energy is E/N = —2-K^br. Finally, the momentum is P/N = (^br, with (' = 27r£, which 
of course is zero for the ground state. 

In the case of the absolute ground state, when there are equal numbers of each type of 
particle so rij = Nj/N = 1/ /, then all limits bj = oo, so b = 11, and the equation can be 
solved by Fourier transforms. The eigenvalues of the integral operators K n are given by the 
Fourier transform of the kernel, 



/oo 
e- lsa K n (a)da = e" |s|/n . (26) 
-oo 



We see K 1 = (K 2 ) 2 - We now define the resolvent operator J by (U + J)(U + K) — (U + 
K)(H + J) = 1. Then 11 + J is also a symmetric matrix of integral operators with difference 
kernels, whose eigenvalues are given by the Fourier transforms 
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[i + J km = %f . S1 h " h( ; 9 7 /2) ,j>™. (27) 

smh(/s/2) smh(s/2) 

The Fourier transform of the densities are then given as 

sinhW-j)/2) 
^' (S) = sinh( /s /2) ' (28) 

confirming rrij = fj(0) = 1 — j/ f so rij = 1/ /. The ground state energy is calculated by 

Parsival's theorem as 

/oo 
e-l'l/ 2 fi(a)(fc = -2[^(1) - V(V/)]//, (29) 
-oo 

with ^(rc) as the digamma function. 

For the excitations about the absolute ground state, there are / — 1 branches, each of 
which is excited pair-wise, as a particle and a hole. The energy and momentum for these 
branches are given as 

^ E M) = ^ ufn Si ^ j//) , .„„ (30a) 
JV / cosh(27ra//)-cos(7rj//)' 

APj(a) = 2 arctan[cot(7rj/2/) tanh(7ra//)] - tt(1 - j//). (30b) 

Eliminating a between the two equations gives a dispersion relation 

Uj (k) = 27r[cos(7ri// - |fc|) - cos(7ri//)]/sin(7ri//), (31) 

for \k\ < 2nj/ f; it is periodic with period 2nj/f. At k = 0, all branches have a common 
velocity t> = 2n/f. This leads to our final result, that if we deviate slightly from equal filling, 
so the concentrations are rij = 1/ f + 5rij, then to leading order in the <5n's, 

AE/N ^X>5, (32) 
J j= i 

and thus the susceptibility exists and is isotropic. 



VII. RESULTS FOR THE HYPERBOLIC MODEL 



We shall make extensive use of the solution to the nearest-neighbor problem to construct 
a solution to the general hyperbolic problem, so let us rewrite the coupled equations for 
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the general problem as £ = (1 + K)bf. The inhomogeneous term is given by Eq. (|i~8l), and 
can be rewritten as £ = £p, by defining £ as an a column vector of integral operators, with 
£j = ^ij^2- This is a natural extension of the previous notation. Since the equations are 
linear, the solution to the general problem is given by superposition as f = rp, r being the 
solution to the nearest-neighbor problem. For the lattice with equal filling, when b = 11, r 
is just a vector of translationally invariant operators. 

Examining the Eqs. (|18|) , ([19]) and (|20|) giving the energy and momentum, we see that 
they can be rewritten as 

Piatt/ N = pltfbr = p ] Q ] br Pl (33a) 
Ei a tt/N = Jtfbr = Jgbrp, (33b) 

and so both can be expressed as quadratic forms in p and e. Here (' = 27r£. The concen- 
trations rij remain unchanged, due to the normalization of p. We see that singularities in 
the groundstate only occur at surfaces of equal filling, the same as for the nearest-neighbor 
case. Again, taking the absolute ground state with equal filling, we can explicitly evaluate 



these expressions using the results of section [VI], giving for the ground state energy 
E latt /N = - f A e(x)dx f A p(x')dx'Re^(l/f + i(x - x')/f) - ^(1 + i{x - x')//)]/tt/. 

J— A J-A 

(34) 

The ground state momentum is of course zero. In Fig. (f|) we show the ground state energy 
as a function of the number of species of fermions /, for several typical densities d, and 
hence various ranges for the hyperbolic interaction. 

For the low-lying excitations, one calculates the shift in the asymptotic momenta, 
weighted by the ground state distributions, and this quantity does not depend upon the 
momentum function p(a). Thus, the only change for the general case comes from the ex- 
pressions Eq. (|33|) for the energy and momentum, so we can write for the excitations 

AP latt = / APj(a-x)p(x)dx, (35a) 
J-A 
f A 

- AEj(a — x)e(x)dx/27r, (35b) 

J-A 
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AE, 



latt 



where APj(a) and AEj(a) are the nearest-neighbor expressions of Eq. fl30|). The corre- 
sponding hydrodynamic velocities for these gapless excitations are all identical, and equal 
to 

v = J-A ex P( 27r x/f)\e(x)\dx ^ 
f I-a exp(27rx/ f)p(x)dx 

This expression is very much like that for a recently solved model for a magnetic fluid, which 
has the Heisenberg-Ising model imbedded in a continuum system ||. In Fig. (||) we show 
the hydrodynamic velocity v as a function of the number of species of fermions /, for several 
typical densities d, and hence various ranges for the hyperbolic interaction. 
The susceptibility can be written using the expression Eq. (|36D for v as 



ae/jv«^£;^ = £x>3. (37) 

Z 3=1 Z 3=1 



Finally, by a thermodynamic argument [HJ, if we introduce fluxes $j = N<pj conjugate to 
the charge of type j, then 

AJB/JV«f E^-?E^?> (38) 
so the charge stiffness D obeys fx -1 = v 2 . 
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FIGURES 

FIG. 1. The dispersion curve e(p) is shown for various values of the density. The en- 
ergy is measured in terms of the nearest-neighbor interaction energy sinh~ 2 (l/(f). In all our 
figures, we have taken as examples, limits of A = 1/2,1,2,4 corresponding to densities of 
d = 0.6735,1.0388,1.7293,3.0596 respectively. In addition, we have the limit A -> 0, d -> 0, 
corresponding to the nearest-neighbor interaction, and the limit A — ► oo, d — > oo, corresponding 
to the 1 /r 2 interaction. 

FIG. 2. The classical ground state distribution function p(x) is shown for representative values 
of the parameter A. 

FIG. 3. The classical ground state distribution function e(x) is shown for representative values 
of the parameter A. It is measured in terms of the nearest-neighbor interaction energy sinh _2 (l/<i). 

FIG. 4. The ground state energy Ei att /N is shown as a function of the number of species of 
fermions /, for representative values of the parameter A. It is measured in terms of the near- 
est-neighbor interaction energy sinh~ 2 (l/d). 

FIG. 5. The hydrodynamic velocity v is shown as a function of the number of species of fermions 
/, for representative values of the parameter A. It is measured in terms of the nearest-neighbor 
interaction energy sinh _2 (l/d). 
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